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ABSTRACT 



To generate the standard microlensing light curve one assumes that the rel- 
ative motion of the source, the lens, and the observer is linear. In reality, 
the relative motion is likely to be more complicated due to accelerations of 
the observer, the lens and the source. The simplest approximation beyond 
the linear-motion assumption is to add a constant acceleration. Microlensing 
light curves due to accelerations can be symmetric or asymmetric depend- 
ing on the angle between the acceleration and the velocity. We show that 
it is possible that some of the previously reported shorter marginal parallax 
events can be reproduced with constant-acceleration models, while the longer, 
multi-year parallax events are ill- fitted by such models. We find that there is 
a generic degeneracy inherent in constant-acceleration microlensing models. 
We also find that there is an equivalent degeneracy in parallax models, which 
manifests itself in short-duration events. The importance of this new paral- 
lax degeneracy is illustrated with an example, using one of these marginal 
parallax events. Our new analysis suggests that another of these previously 
suspected parallax candidate events may be exhibiting some weak binary- 
source signatures. If this turns out to be true, spectroscopic observations of 
the source could determine some parameters in the model and may also con- 
strain or even determine the lens mass. We also point out that symmetric light 
curves with constant accelerations can mimic blended light curves, producing 
misleading Einstein-radius crossing time-scales when fitted by the standard 
'blended' microlensing model; this may have some effect on the estimation of 
optical depth. 



© 2002 RAS 



2 Smith, Mao & Paczynski 

Key words: gravitational lcnsing - Galaxy: bulge - Galaxy: centre - binaries: 
general 

1 INTRODUCTION 

A decade ago Gould (1992) suggested that the Earth's orbital motion may affect long dura- 
tion events of gravitational microlensing, producing a periodic modulation of the event light 
curve. The first such example of this behaviour was detected by Alcock et al. (1995), and a 
number of additional events have been reported by several groups (e.g., Mao 1999; Smith, 
Mao & Wozniak 2002a; Bennett et al. 2002a; Bond et al. 2001). However, most of these 
events were of relatively short duration and showed a strong asymmetry in the light curves, 
rather than a yearly modulation. Two events were long enough to make a yearly modulation 
clearly noticeable: OGLE-1999-BUL-32 = MACHO-99-BLG-22 (Mao et al. 2002; Bennett 
et al. 2002b) and OGLE-1999-BUL-19 (Smith et al. 2002b). 

While the Earth's motion may generate a photometric parallax effect, a similar anomaly 
may be due to binary motion of the source (sometimes referred to as 'xallarap'; Griest & 
Hu 1992; Han & Gould 1997; Paczynski 1997) and/or the lens. In fact such a binary-source 
modulation was detected in the microlensing event MACHO-SMC-1 (Alcock et al. 1997; 
Palanque-Delabrouille et al. 1998; Udalski et al. 1997) and MACHO 96-LMC-2 (Alcock et 
al. 2001). While in these two cases the binary source period was much shorter than the event 
time scale, an alternative situation may also be encountered. In fact, Smith et al. (2002b) 
pointed out that even in the case of a very long event exhibiting an annual modulation of 
its light curve, a model based on a binary source moving with a one year orbital period 
and suitable inclination and orientation of the orbit may generate a light curve identical to 
that expected from the parallax effect. While such a coincidence is very unlikely, it is also 
verifiable in that a binary source should display periodic variations in its radial velocities. 

As most parallax microlensing events are not very long it is likely that the effects they 
show may be well described with a constant acceleration of the Earth, the lens and/or the 
source, with no need to invoke a specific orbit. The aim of this paper is to investigate this 
possibility and to present several consequences of the proposed approach. In Section 2 we 
describe a simple model that incorporates a constant-acceleration term into the microlensing 

* e-mail: (msmith, smao)@jb. man.ac.uk, bp@astro.princeton.edu 
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formalism. We then test the efficacy of this model in Section 3 by fitting a number of 
microlensing events that have been previously identified as exhibiting variations from the 
standard constant-velocity model of microlensing. The events that we have selected were 
all identified as potential parallax microlensing events, and in this section we aim to test 
whether this constant-acceleration model is able to successfully reproduce this behaviour. In 
this section we also mention that there is a generic degeneracy with the constant-acceleration 
model, and that there may be an equivalent degeneracy for parallax models. In Section 4 we 
summarise our work and discuss some of the implications of our findings. In particular, we 
note that the effect of constant-acceleration perturbations is likely to be of little importance 
for binary lenses, and such perturbations are therefore most-likely to be due to accelerations 
of the observer and/or the source. 

2 A SIMPLE MODEL 

In the simplest case of a microlensing event, the light curve is calculated assuming that the 
relative motion of the source, the lens, and the observer is linear and constant (Paczyhski 
1986). The next step in allowing for a more complicated relative motion is to expand it in a 
power series, and to retain the first two terms, i.e. the velocity, v, and the acceleration, a. For 
the purposes of this paper, we define the velocity, v, to be the velocity of the lens relative 
to the observer-source line-of-sight. This acceleration is therefore a relative acceleration, 
but since the following analysis is likely to be valid only for lenses in inertial motion (see 
Section 4), this acceleration can be thought of as the acceleration of the observer and/or 
the source. An acceleration related velocity vector may be defined as v a = atg, where t-& 
is the event time scale: £e = r^/v, te is the Einstein radius and v is the velocity at the 
minimum separation (peak magnification). We expect that any 'anomaly' in the light curve 
will depend on two dimensionless numbers: the ratio of the two velocities: A = v a /v, and 
the angle between them 9. Notice A is simply the physical acceleration in units of r E /t|. 

In general, an acceleration will produce asymmetry in the microlensing light curve, and 
this is the effect which is used to select parallax events. However, if the acceleration is mostly 
in the direction perpendicular to the velocity vector v then there will be no asymmetry, just 
a change in the shape that may be difficult to notice. 

Fig. 1 shows four example light curves with different accelerations. All four cases have 
a minimum impact parameter of u = 0.1, and the source is assumed to move along the 
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x-axis in the absence of acceleration. The dimensionless acceleration parameters along the 
x and y directions, A x and A y , are indicated in each panel. The top left panel shows a 
light curve in which the acceleration is 45° from the x-axis while the bottom left panel 
shows an example where the acceleration is in the same direction as the source motion. As 
expected, both light curves show clear asymmetries. The two right panels show light curves 
in which the acceleration is perpendicular to the rr-axis. The light curves are symmetric in 
such cases. The effect of acceleration may not be noticed, and a standard fit forced upon 
the event may generate incorrect values of the impact parameter and the event time scale; 
we return to this issue briefly in the discussion. A thorough analysis of all such possibilities 
in current experiments is beyond the scope of this paper. Instead, we concentrate on the 
readily available candidate parallax events, in order to check how many of them can be well 
fitted by adding an acceleration vector to the standard model. 

3 ANALYSIS OF CANDIDATE PARALLAX MICROLENSING EVENTS 
USING A CONSTANT ACCELERATION MODEL 

In this section we implement the above constant-acceleration microlensing model and test its 
efficacy by fitting a number of microlensing events. By analysing previously identified candi- 
date parallax microlensing events, we aim to investigate whether this constant-acceleration 
model is able to suitably reproduce the observed (i.e. seemingly parallactic) behaviour. 

3.1 Constant Acceleration Model 

We shall first outline the standard constant-velocity microlensing model. For the standard 
microlensing light curve, the source position (in units of Einstein radius) is given by 

x = t —^- = x (t), y = u , (1) 

where to is the time of maximum magnification and Uq is the corresponding (minimum) 
impact parameter. The magnification is given by (Paczynski 1986) 



where Xo(t) is defined in eq. (1), «A is the physical acceleration in units of r^jt\, and the 
magnification is again calculated using eq. (2). Notice that, since the transverse velocity 
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With acceleration, the (x, y) coordinates must be modified to 

1 1 

x = x (t) + - Acos9 x (t), y = u + - A sin 9 x 2 (t), 



(3) 
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changes, t-& should be understood as the Einstein radius crossing time corresponding to 
the transverse velocity v at time to- Also, in general, uo may no longer rigorously be the 
minimum impact parameter, but may still provide a useful approximation, particularly if A 
is small. 

To fit a microlensing light curve with the standard model, one needs a minimum of four 
parameters, uo, to, £e, and the baseline flux (or magnitude), fo- One also often requires an 
extra blending parameter, which accounts for the additional light contributed by the lens 
and/or other stars in the crowded field. To do this, we use a parameter, / s , defined as the 
fraction of light contributed by the lensed source at the baseline. To incorporate the constant 
acceleration, we need two extra parameters, A and the angle 9 (or equivalently, A x and Ay). 

We look for any degeneracy in the parameter space for this constant-acceleration model 
by fitting each event 1000 times, with the initial parameter guesses selected from wide range 
of randomly chosen values. We find that two sets of parameters exist for each \ 2 minimum. 
The parameters are identical except for the time-scale, £e, and the acceleration parameters, 
A and 9. However, the magnitude of the physical acceleration in units of rs/day 2 (i.e., 
|a/r E | = Ajt\, as defined in §2) is the same for both fits, although the direction of this 
vector (i.e., 9) differs. In Appendix Al, we show that this degeneracy can be understood 
analytically. An example of this degeneracy is illustrated in the inset of Fig. 2 for the 
microlensing event sc6_2563; even though the two trajectories look very different they both 
result in identical light curves. The full analysis of the event sc6_2563 is described in Section 
3.2.1. 

It should be noted that one may expect this degeneracy to manifest itself in parallax 
models. In particular, this could be important for short- duration events that exhibit weak 
parallax signatures, since in this case the Earth's acceleration can be approximated to be 
constant (i.e. the regime in which this degeneracy is expected to occur). By testing this 
hypothesis on the event sc6_2563, we are able to verify the existence of this degeneracy. 
However, it seems that this may only become apparent for weak parallax candidates, since 
a more convincing parallax candidate (sc33_4505; Section 3.2.3) shows no such degeneracy 
(see Appendix A). A full discussion of this degeneracy, including the application to events 
sc6_2563 and sc33_4505, can be found in Appendix A2. 
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Table 1. The best standard model (first row), the best parallax models (second and third rows), and the best constant- 
acceleration model (fourth and fifth rows) for the OGLE-II event sc6_2563. The final column shows the \ 2 an d the number of 
degrees of freedom (dof) for each model. The parameters are explained in §2 & §3. 8 is given in units of radians. There are two 
sets of constant-acceleration parameters since each \ 2 minimum has two degenerate fits with differing values for (e, A and 9. 
There are also two sets of parallax parameters, owing to the parallax degeneracy that is analogous to this constant-acceleration 
degeneracy. Both degeneracies are discussed in Section 3.1 and Appendix A. 
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3.2 Analysis of candidate parallax microlensing events 

We proceed by analysing a selection of events that have previously been identified as can- 
didate parallax microlensing events. The parallax model used here requires 7 parameters: 
the 5 parameters from the standard model plus two additional parameters to describe the 
lens trajectory in the ecliptic plane, namely the Einstein radius projected onto the observer 
plane, -te, and an angle ip in the ecliptic plane, ip is defined as the angle between the helio- 
centric ecliptic rr-axis and the normal to the trajectory (This geometry is illustrated in Fig. 5 
of Soszyhski et al. 2001). These candidate parallax events have been selected to provide a 
broad range of parallactic signatures, ranging from marginal cases to those displaying more 
prominent effects. 



3.2.1 sc6_2563 

This event was detected by the OGLE collaboration (Udalski et al. 2000). Its asymmetric 
nature was identified during a parallax search of the OGLE-II database (Wozniak et al. 
2001), reported in Smith et al. (2002a). It was classified as a 'marginal' candidate, since the 
deviations from the standard model are not particularly pronounced and the improvement 
in x 2 is on ly slight. The best-fit parameters for the standard and parallax models are given 
in Table 1, and the corresponding light curves are plotted in Fig. 2. The inset of this figure 
shows the two degenerate trajectories for the constant-acceleration model; as was discussed 
above, each x 2 minimum has two degenerate fits. Both of these trajectories produce iden- 
tical light curves, even though the parameter values for £e, A and 9 differ. This table also 
includes an additional parallax fit, which has a x 2 difference of 0.8 compared to the best-fit 
parallax model. This is a manifestation of the parallax degeneracy that is analogous to this 
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constant-acceleration degeneracy. See Section 3.1 and Appendix A for a discussion of these 
degeneracies. 

It would be expected that a constant-acceleration fit may be suitable for this event 
since the duration is not very long, and the asymmetry shows no signs of modulation. 
Indeed, when this event is fit with the above constant-acceleration model the best-fit x 2 value 
improves upon the best-fit parallax value from 347.4 to 332.6. If the error bars are rescaled 
so that the best-fit x 2 P er degree of freedom is renormalised to unity, then the difference 
in x 2 between these two fits is 9.5 for no additional free parameters, i.e. a significant 3-cr 
improvement. The best-fit parameters and corresponding light curve for the acceleration 
model are given in Table 1 and Fig. 2, respectively. From this figure it is clear that the 
parallax and constant-acceleration models produce similar fits. However, both models seem 
to over-predict the flux around 1000 < t < 1150 days, implying that there may be some 
non-constant contribution to the acceleration that is not parallactic in nature. Notice that 
the standard microlensing model predicts a blending parameter / s m 0.8 while other models 
predict / s ~ 1. Correspondingly, the Einstein radius crossing time is about 25% larger for 
the standard model than for the other models; this has implications for the optical depth 
estimate in the experiments (see Section 4). 

3.2.2 OGLE-1999-CAR-1 

This event was detected in real-time by the OGLE Early Warning Systemt (Udalski et al. 
1994) toward the Carina spiral arm. It has been found to exhibit systematic deviations from 
the standard model (Mao 1999). The duration for this event is well over 100 days, and it 
was concluded that these deviations are due to the parallax effect. 

As with sc6_2563, the deviations from the standard model for OGLE-1999-CAR-1 show 
no signs of modulation, and so one would expect that this event may be suitably approx- 
imated by a constant acceleration model. This appears to be correct, since the best-fit 
constant-acceleration model is able to prove a very similar fit (with a slightly worse x 2 value 
of 619.4 compared with the parallax best-fit value of 614.8). The best-fit parameters for all 
models are given in Table 2, and the corresponding light curves are shown in Fig. 3. It is 
interesting to note that both of the degenerate acceleration fits have values of A = v a /v 
greater than 1 (A — 1.4, and A = 2.4). This implies that v a , the change in velocity due to 

t http://www.astrouw.edu.pl/~ogle/ogle3/cws/cws.html 
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Table 2. The best standard model (first row), the best parallax model (second row), and the best constant-acceleration model 
(third and fourth rows) for OGLE-1999-CAR-1. The parameters arc explained in §2 & §3. 8 is given in units of radians. There 
are two sets of constant-acceleration parameters since each \ 2 minimum has two degenerate fits with differing values for <b, A. 
and 8 (see Section 3.1 and Appendix Al). 
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this constant acceleration in a time £e, is greater than the v, the velocity at the point of 
closest approach. 

Also of note is the difference in t E among the three models. The values range from t E = 
126 days for the parallax model, to as much as t E = 190 days for the constant-acceleration 
model. These differences in the time-scale are closely connected with their differences of the 
blending parameter, due to potential degeneracy between these two parameters (Wozniak 
& Paczynski 1997). 



3.2.3 sc334505 

This OGLE-II event was also found during the search of Smith et al. (2002a), but was 
classified as a 'convincing' parallax candidate. Unlike sc6_2563, the duration of this event 
was much longer (t E ~ 200 days, cf. t E ~ 100 days for sc6_2563), and the data show clear 
signs of deviation from the standard model in both the rising and declining branch of the 
light curve. Therefore, if this asymmetry is indeed due to the parallax effect then one would 
suspect that a constant-acceleration model would be unable to reproduce this behaviour 
(since the acceleration of the Earth cannot be considered constant during this time span). 
On the other hand, if this asymmetry is due to other causes, the constant-acceleration model 
may provide a better fit. 

The best-fit parameters for the standard and parallax models are given in Table 3, and 
the corresponding light curves are shown in Fig. 4. A problem arises when this event is fit 
with the constant-acceleration model, in that two minima are identified, with x 2 — 187.3 
and x 2 — 227.3, respectively. Due to the degeneracy discussed in §3 and Appendix Al, there 
are two degenerate fits corresponding to each of these two x 2 values. The two degenerate 
fits that correspond to x 2 = 187.3 have extremely large time-scales (t E = 1812.8 day and 

= 4880.5 day), and the lensed star contributes only 2.3 per cent of the total baseline 
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light (i.e., f s = 0.023). These values do not seem physical to us, and therefore the fits 
corresponding to this x 2 minimum are disregarded. The best-fit parameters for the two 
degenerate fits corresponding to x 2 — 227.3 are more 'feasible'; they are given in Table 3 
and the corresponding light curve is shown in Fig. 4. 

However, the discarded 'unfeasible' constant-acceleration models have an overall best x 2 
value (187.3). Since this x 2 value is better than the best-fit parallax model, it highlights 
that there could be a deficiency in the parallax fit. This deficiency in the parallax fit can be 
seen in the residual plot in Fig. 4; systematic residuals from the parallax model can clearly 
be seen, and the remaining 'feasible' constant-acceleration fit is also unable to successfully 
reproduce this behaviour. 

These residuals from the parallax model could imply that this event is exhibiting some 
weak binary-source signatures (see, for example, Alcock et al. 2001). We performed a pre- 
liminary investigation to test this hypothesis, fitting the event with a simple binary-source 
model. The basic details of this model can be found in Smith et al. (2002b), but the version 
used here has two differences from the earlier work: firstly, it incorporates the parallax mo- 
tion of the Earth (which, in most cases, should have an effect on the light curve since the 
duration of this event is well over one yearf); and secondly, only face-on elliptical orbits are 
considered. 

The best binary-source plus parallax fit that we are able to identify has a much better 
X 2 value (x 2 = 162.4) compared to all of the previous models. However, we find a range of 
possible binary-source fits with x 2 l ess than 180. The single best-fit solution is presented in 
Table 4 and Fig. 4. This model predicts a very large value for the Einstein radius projected 
into the observer plane (f E = 108 au), which implies that the parallax motion of the Earth 
has little effect on the light curve, i.e., the majority of the perturbations from the standard 
model arise due to the binary nature of the source. It seems that this model is able to fit 
the systematic residuals from the parallax model, although it predicts peculiar bumps in 
the light curve (most notably during the break between observing seasons). We present this 
model to show that such fits are possible, but we do not attempt a more detailed analysis 
of this event in this work. 

If this event is indeed affected by the binary motion of the source then a spectroscopic 

t It should be noted that the parallax effect is expected to be negligible in the case of binary-source microlensing when the 
lens and sources both lie in the bulge, so-called bulge-bulge lensing. 
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Table 3. The best standard model (first row), the best parallax model (second row), and the best constant-acceleration model 
(third and fourth rows) for the OGLE-II event sc33_4505. The parameters are explained in §2 & §3. 9 is given in units of 
radians. There are two sets of constant-acceleration parameters since each \ 2 minimum has two degenerate fits with differing 
values for 4e, A. and 9 (see Section 3.1 and Appendix Al). 

Model t t E (day) ii„ /„ /, V ■* 6 X 2 /dof 

S 647.22 ± 0.32 165.9+? ^ 0.412 + 55- iSo 1146.3 ± 1.5 0.637+9'?™ 490.1/179 

217.4/177 

642.36 ± 0.46 115.38+^ ** 0.5596 ± .0014 1149. 4 + J ^ 1.000+° Q&i — °- 561 io o24 4 252 + 0.015 227.3/177 

— 1 ' 174 lo o93 2 268 ± 029 227.3/177 



Table 4. The best-fit parameters for the elliptical binary-source model for the OGLE-II event sc33_4505. The first seven 
parameters correspond to the usual parallax microlensing parameters (see §3), and the final six describe the binary nature of 
the source. The mass (flux) ratio, M. {T), is defined as the mass (flux) of the first binary source divided by the total mass 
(flux) of the two sources. A detailed description of the binary-source model can be found in Smith ct al. (2002b). The errors 
have been omitted since they were found to be misleading, owing to the complexity of the x 2 surface. 

*0 *E (day) u fo fs f-E 6 

675.2 140.5 0.573 1143.2 1.00 107.9 5.623 

orbital period (day) semi-major axis (r^) M T eccentricity phase x 2 /dof 

159.48 1.057 0.105 0.192 0.943 1.216 162.415/71 

followup should find the source to be a spectroscopic binary, and this would reduce the 
number of free parameters in our model. Furthermore, if the projected separation of the 
two sources can be estimated, then this would enable us to make a direct determination of 
the lens mass (Han & Gould 1997). The source baseline magnitude (including blending) is 
about / = 15.7 mag, and so it should be within easy reach of modern large telescopes. 

3.2.4 OGLE-1999-BUL-19 and OGLE-1999-BUL-32/MACHO-99-BUL-22 

These two events display unambiguous deviations from the standard model. They were de- 
tected by the OGLE Early Warning System (with OGLE-1999-BUL-32 being independently 
detected by the MACHO collaboration's Alert System§, by which it was named MACHO- 
99-BUL-22). Both of these events also appear in the difference image analysis catalogue of 
Wozniak et al. (2001), and are labelled sc40_2895 and sc33_3764, respectively. For the pur- 
poses of this paper, we fit these events using the OGLE difference image analysis data from 
Wozniak et al. (2001). 

OGLE-1999-BUL-19 (Smith et al. 2002b) displays prominent multiple peaks that are 
well fit by a parallax model, while OGLE-1999-BUL-32 (Mao et al. 2002; Bennett et al. 

§ http: / / darkstar.astro.washington.edu/ 
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Table 5. The best \ 2 values for the two unusual events OGLE-1999-BUL-19 (Smith ct al. 2002b) and OGLE-1999-BUL-32 
(Mao et al. 2002; Bennett et al. 2002b). 

OGLE-1999-BUL-19 OGLE-1999-BUL-32 

Standard x 2 / dof 17910/312 576.3/264 

Parallax x 2 /dof 590.1/310 278.2/261 

Constant Acceleration x 2 /dof 11420/310 464.4/261 

2002b) is the longest duration event ever discovered. Both events are suitably fit by the 
parallax model, and both show significant improvement in x 2 between the standard and 
parallax models. 

Since the microlensing amplification for both of these events lasts well over two years, we 
expect that a constant-acceleration model may not be adequate to describe their behaviour. 
By fitting with the constant-acceleration model we find that this is true; neither of these 
events are well fit. The best-fit x 2 values for these events are given in Table 5. For OGLE- 
1999-BUL-19, the x 2 f° r the constant-acceleration model is 11420, which is far worse than 
the parallax model of 590.1. Similarly, for OGLE-1999-BUL-32, the x 2 f° r the constant 
acceleration model is 464.4 compared with 278.2 for the parallax model. 

4 DISCUSSION 

In this paper, we have illustrated the effect of acceleration on microlensing light curves. We 
studied whether the known parallax microlensing events can be equally fitted by a model 
with a constant acceleration. We found that for the shorter marginal parallax microlensing 
events, the constant acceleration is indeed sufficient. However, for the longest events (notably, 
OGLE-1999-BUL-32 and OGLE-1999-BUL-19), the constant acceleration fails to provide a 
satisfactory model, so their parallax nature seems secure. 

It is important to examine whether there is an empirical way of separating 'parallax 
events' from 'microlensing events with acceleration'. For binary lenses it can be shown that 
such constant-acceleration perturbations are unlikely to be of importance. This can be un- 
derstood by considering three regimes of binary lenses separately: close binaries (i.e., those 
with separation b <C te), intermediate binaries (6 ~ te), and wide binaries {b ^> te). One 
would expect sufficiently large accelerations for close separation binaries. However, unless 
the impact parameter is extremely small (i.e. the source passes very close to the centre of 
mass of the lens), it is unlikely that the resulting light curve will be differentiable from a 
standard single- lens light curve (Gaudi & Gould 1997). Even in the rare cases where such 
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small impact parameters occur, one would expect the primary deviations to be caused by the 
binary structure, rather than the acceleration (for example, the binary-lens event MACHO- 
99-BLG-47, Albrow et al. 2002). For intermediate binaries where the separation is of the 
order of the Einstein radius, one would also expect the primary deviations to arise from the 
binary structure, rather than the acceleration. Even in the absence of caustic crossings, the 
effect of acceleration should be relatively small compared to the distortion of the light curve 
due to the other mass (For example, typical lens geometries - with the source lying in the 
bulge at a distance of approximately 8 kpc and the lens lying half-way in between - predict 
te ~ 4 au (Miens/il^,) 1 / 2 , where Mi ens is the total lens mass; therefore, for face-on orbits with 
separations of approximately 1 r E , a typical orbital period is P ~ 8 yr (Mi ens / 'Mq) 1 / 4 , which 
implies that the resulting acceleration should produce only weak perturbations to the light 
curve compared to the distortion from the companion mass). For the case of wide separation 
binaries, it is clear that the accelerations should be small compared to the Earth's accelera- 
tion, since binaries with separation much greater than r E are expected to have orbital periods 
of longer than 8 years. Therefore, we conclude that the effect of constant accelerations in 
binary lenses is unlikely to be of importance. 

In principle, there is a simple test to determine whether an acceleration is due to binary 
motion of the source or the observer: for the case where the acceleration is induced by the 
source's binary motion, the source must show (periodic) changes in radial velocity, which 
can be verified even after the event is over. In Appendix B, we illustrate the expected radial 
velocities for the simple case of circular binary orbits. We find that the expected radial 
velocity amplitude is of the order of 35kms _1 ^4 1//4 sin?/ cos 1//4 2, where % is the inclination 
angle; there is also some weak dependency on t E , the transverse velocity of the lens, v, and 
the masses. The expected orbital period depends on the unknown lens transverse velocity 
and the inclination angle, but is of the order of ~ years. This radial velocity can, in principle, 
be measured. 

In the right two panels of Fig. 1, we showed that the light curves are symmetric when the 
acceleration is perpendicular to the source trajectory. In such cases, the effect of acceleration 
may be difficult to notice. This may even be true for some cases for which the acceleration is 
not exactly perpendicular to the source motion, particularly since the real light curves have 
significant gaps due to bad weather, etc. If we force the standard fit to the microlensing 
light curve, the time-scale obtained from the fit may be different from the true value. To 
illustrate the effect, we generate artificial light curves using Monte Carlo simulations. We 
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adopt tft = 30 days, and a sampling frequency of one observation per day; the observations 
lasts from — 2£e to 2£e centred around the peak of the light curve. The observational errors 
are assumed to be 0.05 magnitudes at the baseline and scale according to photon Poisson 
noise. We find that a standard fit incorporating blending is significantly better than a fit 
without blending. In other words, a microlensing light curve with constant acceleration can 
mimic light curves with blending. Interestingly, the fit with blending under-estimates the 
true t-E by 20% for both acceleration values (A = ±1.5) shown in the right panels of Fig. 
1. This can have important implications for the calculation of optical depth. The optical 
depth, r, is estimated from experiments using the following formula (see, for example, Han 
& Gould 1995), 
N + 

r = T ' l (4) 

2W.r,^e(* E) i)' yJ 

where iV* is the total number of monitored stars, T* is the experiment duration in years, 
e(tE,i) is the detection efficiency, and t^ i is the Einstein radius crossing time for the i-th 
event (i = 1, . . . ,N). Obviously, if the time-scale is underestimated then this will lead to 
an underestimation in the optical depth, and vice versa. However, a quantitative analysis 
requires a detailed knowledge about the binary parameters of the lens and the source, and 
sampling frequencies in observations. This is clearly an area that deserves further study. 
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APPENDIX A: DEGENERACIES 

In this appendix we first describe the degeneracy that is found to arise in the constant- 
acceleration model. This is followed by a discussion of an equivalent degeneracy in the 
parallax model, which can occur for short duration events where the Earth's acceleration 
may be approximated as constant. 

Al Degeneracy in Constant- Acceleration Models 

In Section 3 we discuss the degeneracy inherent in fitting microlensing light curves with 
the constant-acceleration model. It is found that each x 2 minimum has two corresponding 
fits, each with differing values of time-scale, t^, and the acceleration parameters, A and 9. 
However, the magnitude of the physical acceleration in units of rs/day 2 (i.e., |a/rE| = A/t^) 
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is the same for both fits, although the direction of this vector (i.e., 6) differs. In this subsection 
we show analytically why this degeneracy occurs. 

Given a combination of parameters (£e, uo, etc), the magnitude of the source displacement 
in the lens plane (in units of the Einstein radius) at a time (t + St) can be determined from 
equations (1-3), 

u 2 (t + St) = (^a cos 9 + + (^a sin 9 ^y- + u ^ , (Al) 

where a = |a/r E | = A/t\\ is the magnitude of the physical acceleration in units of the 
Einstein radius. 

This equation can be expanded to give, 

u 2 (t + St) = a 2 ^- 

f^cos9\ o 

+ a — )<*> 



+ ^^—2 + -u asin#^ (Stf 



J E 2 

+ «g. (A2) 

For two degenerate solutions i and j, the magnitude of the source displacement, Ui(to + St) 
and Uj(to + St), must be the same for all values of St (since this is required to produce 
identical light curves), i.e., u 2 (t + St) = u 2 (t + St) for all St. Since this is true for all 
values of St, using equation (A2) we can equate powers of St. This gives the following four 
equations, 

hi = a,j (A3) 

„ cos 9i „ cos 9j , . 4 . 

a,—^ = aj — (A4) 



1 „ . „ 1 

2 

] ,3 



+ u^hi sin 9i = + u 0J aj sin 9j (A5) 



«o,i = <r ( A6 ) 
From equations (A3) and (A6) we can see that both a and uq must be constant for degenerate 
solutions (the sign of u is irrelevant due to the symmetry of the acceleration model; by 
taking only u > one still obtains the full set of solutions). Also, from equations (A4) and 
(A5) we find that the quantities [cos#/£ E ] and [1 / (t E ) 2 + u a sin 9] must also be constant for 
degenerate solutions. Therefore, the following two equations must hold for every degenerate 
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solution, 

-r- = ci, (A7) 
— + m a sin = C2, (A8) 

^E 

where Ci and c 2 are constants defined for any given x 2 minimum. The parameter 9 can then 
be eliminated from these two equations to give the following expression for £e, 

f(w) = (auodfw 3 + (cj - a 2 ul)w 2 - 2c 2 w + 1 = 0, (A9) 

where w = t\. We see that the above equation is a cubic equation in w. As the coefficients 
are all real, the number of real solutions is either one or three. Recall that since w = t|, 
physical solutions to equation (A9) must clearly correspond to positive values of w. 

As /(— oo) = — oo, and f(w = 0) = 1, we must have one negative real solution for 
w, which is unphysical and should be discarded^. It is easy to show that the number of 
positive (i.e. physical) solutions for w must be either zero or twoll. This follows because 
f(w = 0) = 1 and f(w = +oo) = +oo. If f(w) is everywhere positive for w > 0, then 
the number of positive solutions is obviously zero. If, however, f(w) is negative at some 
point Wi(> 0), then there must be one positive solution in the region to Wi, and another 
positive solution in the region w\ to +oo. This therefore immediately implies that if there 
is one good constant-acceleration fit (i.e. one physical solution for w), then there must be 
another degenerate fit with the same x 2 , just as we have found numerically. These two fits 
have identical values of u and a = A/t^, but different values of A, 9 and t^- 

A2 A new parallax degeneracy 

A consequence of this constant-acceleration degeneracy is that a similar degeneracy should 
arise for short- duration events that exhibit weak parallax signatures. This is because the 
acceleration in the parallax model can be approximated as constant provided that the du- 
ration of the event is sufficiently short. In the following subsection we consider this new 
degeneracy and provide an example using the event sc6_2563, which is introduced in Section 
3.2.1. 

There are currently well-known degeneracies inherent in parallax microlensing: for ex- 

» It is possible to have three negative solutions to this equation, although this possibility can be ignored since this would result 
in no positive (i.e. physical) solutions 

II Although in rare cases, the two positive solutions can be identical 
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ample, Gould, Miralda-Escude & Bahcall (1994) discussed the continuous degeneracy that 
occurs for events with extremely weak parallax signatures where only one component of the 
relative lens velocity, v, is measurable; another four-fold discrete degeneracy was identified 
in the case of satellite parallax measurements (Refsdal 1966; Gould 1994), although it was 
shown that this problem can be resolved using higher-order effects (Gould 1995). 

As we have discussed above, in the case of the weak candidate parallax events, such as 
sc6_2563 (see section 3.2.1), the best constant-acceleration fit mimics the best parallax fit. 
We also found that each x 2 minima for the constant-acceleration model possesses two degen- 
erate fits, each with identical values for impact parameter, uo and acceleration magnitude, 
a = A/t\, but different values of A, 9 and t E . However, since the above constant-acceleration 
formalism is rotationally invariant, it is possible to rotate one of these degenerate constant- 
acceleration fits so that both degenerate fits have their acceleration vectors co-aligned. There- 
fore, if one of these best-fit constant-acceleration solutions mimics the existing parallax fit, 
then the degenerate counterpart (once rotated through the required angle so that the accel- 
eration is also in the same direction as the Earth's acceleration) should likewise mimic an 
equivalent parallax fit. Since this new parallax fit must be different from the existing one, 
this implies that there should be two degenerate parallax fits. Obviously, this argument only 
holds in the regime of constant-acceleration, i.e., this would not apply to parallax events 
where the duration is sufficiently long to invalidate the constant-acceleration approximation. 

We illustrate this new degeneracy by providing an example using the event sc6_2563, 
which was introduced in Section 3.2.1. To identify any degenerate solutions in the parallax 
parameter space we fit this event 1000 times, taking random values for the initial parameter 
guesses. From this we identify another solution in the parameter space that has a slightly 
worse x 2 value than the previously identified best-fit parallax solution. Both fits are presented 
in Table 1. If the x 2 values are renormalised to enforce the x 2 P er degree of freedom to 
be unity for the best-fit parallax model, then the difference in x 2 is l ess than 0.5. Unlike 
the two degenerate constant-acceleration fits, these parallax fits have different values for 
the magnitude of the impact parameter, uq. However, this is due to the difference in the 
definition of uq in the two models. 

Figure 5 shows the trajectories in the lens plane for the two parallax fits, along with the 
two degenerate constant-acceleration fits. Note that the A' trajectory has been rotated so 
that its acceleration is pointing in the same direction as the A trajectory (the rotation does 
not affect the light curve). Although each parallax trajectory deviates from its corresponding 
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constant-acceleration counterpart, it can be seen that the distance from the origin (i.e., the 
quantity that determines the magnification; the parameter u in eq. 2) is very similar. As 
expected, this figure shows that the difference between the constant-acceleration and parallax 
trajectories only becomes apparent away from the point of closest approach, i.e. when the 
direction of the Earth's acceleration has changed. This corresponds to the wings of the light 
curve, but since the magnification is only weak in this region the resulting deviations in the 
light curve are difficult to detect. 

We also investigated the event sc33_4505, which is introduced in Section 3.2.3. This event 
exhibits more prominent deviations, and the constant-acceleration model is found to provide 
a significantly worse fit than the parallax model (see Table 3). Again, we fit this event 1000 
times to try and identify any degeneracy that might occur for the parallax model, but we 
are unable to find any such degenerate fits. The next-best parallax fit has x 2 — 227.9, which 
is more than 3 a away from the best-fit value of \ 2 = 217.4. Therefore, we conclude that 
there appears to be no such parallax degeneracy for this event, owing to the fact that the 
deviations in the light curve are relatively prominent, i.e., these deviations are unable to be 
reproduced by a constant-acceleration fit. 

APPENDIX B: RADIAL VELOCITY 

Here we estimate the expected radial velocity for the case where the acceleration is induced 
by the binary motion of the source. For simplicity, we consider the case where both orbits 
are circular. Let us denote the mass of the lensed source star as M and that of its binary 
companion as M', the separation between the stars as d, and the inclination of the orbital 
plane as i (i — 90° implies an edge-on orbit). 

In the centre-of-mass rest frame, the maximum radial velocity is given by 



The physical acceleration of the lensed source due to its companion is given by GM'/d 2 . 
The component perpendicular to the line-of-sight is a = (GM'/d 2 )(l — sin 2 sin 2 i) 1 ^ 2 ; here 
4> is the angle of the position of the lensed star from the line where the orbital plane and 
the plane of sky intersect. Notice that our constant acceleration assumption requires that 
does not change significantly during the lensing event, i.e. for circular orbits A0 << 1 




(Bl) 
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implies that the period of the binary P » 2-Kt^. Combining the expression of a with the 
definition of A, we have 

Cancelling out d from eqs. (Bl) and (B2), we obtain the maximum radial velocity 
(M'GvV 1 * k 1/4 I M' V' 2 sim 

-=(— ) A '{mTw) (l _ 3illWtr (B3 » 

where v(— r-^/t^) is the lens transverse velocity perpendicular to the observer-source line. 

Numerically, we have 

X (W lOOdays v \ 1/4 , 1/4 / W \ 1/2 sini 

w r = 35kms _1 — A 1/4 77^-(B4) 

\M Q t E lOOkms-V {m + M'J (i_ sin 2 0sin 2^V^ ] 

Notice that the dependence on the masses, A and the lens transverse velocity, v, are fairly 
weak. 

From Kepler's third-law, the period of the binary is given by P = 27rd 3 / 2 G~ 1 / 2 (M + 
M') -1 / 2 . Substituting d from eq. (B2) into P, we obtain 



r / M \ / *e lOOkms- 1 ^ 374 , V4 M> \ 1/2 / . 2j . 2 .\ 3 / 8 

P = 221 days ^ A~ 3/4 ( 1 - sin 2 0sin 2 i) 

y \M & ) \100days v J \M + M' J V y J 

The dependence of period on the unknown lens transverse velocity and the inclination is 
fairly sensitive. For events that show constant accelerations, the transverse velocity may be 
fairly low and hence the period can be of the order of years. 
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Figure 1. Four example microlcnsing light curves incorporating acceleration effects. All four panels assume the same minimum 
impact parameter uq = 0.1 in the case of no acceleration. The accelerations along the x and y directions are labelled at the 
top right of each panel. The source is assumed to move along the a;-direction in the case of no acceleration. In all panels, the 
dashed and solid lines show the light curves with and without acceleration, respectively. 
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Figure 2. The /-band light curve for OGLE-II event sc6_2563 from difference image analysis, with flux given in units of 10 
ADU. The data for the first season has been omitted from this plot. The /-band baseline magnitude for this event is 17.7, 
rising to 15.1 at the peak. The top panel shows the residual flux (the observed data points with the parallax model subtracted). 
This figure shows that the constant-acceleration model and the parallax model produce very similar fits, with the former model 
providing a slight improvement in \ 2 - However, both models appear to over-predict the flux around 1000 < t < 1150 days. The 
inset illustrates the two possible trajectories for the constant-acceleration model, described by the two fits A (bottom) and A' 
(top) given in Table 1. The dotted circle represents the size of the event's Einstein radius, rE, and the axes are in units of this 
Einstein radius. Both trajectories result in identical light curves, even though the parameter values for tj;, A and 8 differ (see 
Section 3 and Appendix Al for a discussion of this degeneracy). 
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Figure 3. The /-band light curve for event OGLE-1999-CAR-1 from difference image analysis. The data for the first season 
have been omitted from this plot. The top panel shows the residual magnitude (the observed data points with the parallax 
model subtracted). This shows that the parallax model and the constant-acceleration model produce almost identical fits. 
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Figure 4. The /-band light curve for the OGLE-II event sc33_4505 from difference image analysis. The data for the third season 
have been omitted from this plot. The /-band magnitude for this event is 15.7 at the baseline and rises to 15.0 at the peak. 
The top panel shows the residual flux (the observed data points with the parallax model subtracted). It is clear from this panel 
that there are systematic residuals from the parallax fit. The constant acceleration model appears to fit the data much better, 
but this still seems to be deficient. However, the 'binary-source plus parallax' model provides the best-fit (see Section 3.2.3). 
Spectroscopic observations could confirm whether this event is being affected by binary motion of the source; such observations 
would enable constraints to be put on some of the model parameters, and may also result in a direct measurement of the lens 
mass. 
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Figure 5. The degenerate parallax (solid lines) and constant-acceleration (dashed lines) trajectories for the OGLE-II event 
sc6_2563. These trajectories are described in the lens plane and correspond to the paths taken by the lens relative to the 
observer-source line-of-sight (i.e. the origin). Each degenerate fit is labelled according to Table 1. Note that the A' trajectory 
has been rotated so that its acceleration is pointing in the same direction as the A trajectory (the rotation does not affect the 
light curve). The contours of constant amplification are represented by the dotted circles (since the amplification is dependent 
on the separation of the lens from the observer-source line-of-sight; see eq. 2), and the Einstein radius is given by the solid 
circle. The inset shows the light curve for this event. The crosses on each trajectory denote 25 day intervals, as labelled on the 
inset light curve. Note that all of the trajectories presented here result in (practically) the same light curve. 
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